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We have analyzed the latest NA50 data on J/ip suppression in Pb + Pb collisions. J/ip production 
is assumed to be a two step process, (i) formation of cc pair, which is accurately calculable in QCD 
and (ii) formation of J/ip meson from the cc pair, which can be conveniently parameterized. In 
a pA/AA collision, the as the cc pair pass through the nuclear medium, it gain relative square 
momentum at the rate of e 2 per unit path length. As a result, some of the cc pairs can gain 
enough momentum to cross the threshold to become an open charm meson, leading to suppression 
in pA/AA collisions. The parameters of the model were fixed from experimental data on the total 
J/ip cross section as a function of effective nuclear length. The model without any free parameter, 
give excellent description of NA50 data on Et dependence of J /if) to Drell-Yan ratio. The model 
was applied to predict the Et dependence of J/ip at RHIC energy. Much larger suppression of J/ip, 
O i 1 in agreement with other model calculations are predicted. 
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I. INTRODUCTION 

In relativistic heavy ion collisions J/ip suppression has been recognized as an important tool to identify the possible 
phase transition to quark-gluon plasma. Because of the large mass of the charm quarks, cc pairs are produced on 
a short time scale. Their tight binding also makes them immune to final state interactions. Their evolution probes 
the state of matter in the early stage of the collisions. Matsui and Satz jlj predicted that in presence of quark-gluon 
H plasma (QGP), binding of cc pairs into a J/ip meson will be hindered, leading to the so called J/ip suppression in 
heavy ion collisions [jj]]. Over the years several groups have measured the J/ip yield in heavy ion collisions (for a 
review of the data and the interpretations see Refs. In brief, experimental data do show suppression. However 

O ^ this could be attributed to the conventional nuclear absorption, also present in pA collisions. 

The latest data obtained by the NA50 collaboration [Q on J/ip production in Pb+Pb collisions at 158 A GeV 
is the first indication of the anomalous mechanism of charmonium suppression, which goes beyond the conventional 
suppression in a nuclear environment. The ratio of J/ip yield to that of Drell-Yan pairs decreases faster with Et in the 
• • ■ most central collisions than in the less central ones. It has been suggested that the resulting pattern can be understood 
m ^ ' in a dcconfinement scenario in terms of successive melting of charmonium bound states (J] . In a recent paper , Blaizot 
et al. [jp have shown that the data can be understood as an effect of transverse energy fluctuations in central heavy 
ion collisions. Introducing a factor e = Et/Et(V) , assuming that the suppression is 100% above a threshold density 
(a parameter in the model), and smearing the threshold density (at the expense of another parameter) the best fit to 
the data was obtained. Extending the Blaizot's model to include fluctuations in number of NN collisions at a fixed 
impact parameter, NA50 data could be fitted with a single parameter, the threshold density, above which all the J/ip 
mesons melt [||. Assumption that all the J/ip mesons melt above a threshold density, implicitly assume that QGP like 
environment is produced in the collision. NA50 data could be explained in the conventional approach also, without 
invoking QGP like scenario. Capella et al. JjJ analyzed the data in the comover approach. There also, the comover 
density has to be modified by the factor e. Introduction of this adhoc factor e can be justified in a model based on 
excited nucleons represented by strings || . 

Aim of the present paper is to show that while in conventional approach, nuclear suppression is not sufficient to 
explain NA50 data, the data are very well described in a model of Qiu, Vary and Zhang ||, where the suppression 
due to nuclear environment is treated in an unconventional manner. 



II. MODEL 



Recently, Qiu,Vary and Zhang || proposed a model to describe the J/ip suppression in nucleon- nucleus/nucleus- 
nucleus collisions. For the sake of completeness, we will briefly describe the model. Qiu, Vary and Zhang assumed 
that the production of J/ip meson is a two step process, (i) production of cc pairs with relative momentum square q 2 , 
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and (ii) formation of J/tp mesons from the cc pairs. Step (i) can be accurately calculated in QCD. The second step, 
formation of J/tp mesons from initially compact cc pairs is non-perturbative. They used a parametric form for the 
step (ii), formation of J/tp from cc pairs. The J/tp cross section in AB collisions, at center of mass energy y/s was 
then written as, 

°A+B^J/il>+x(s) = j d 1 2 f aab Q2 ^ J J dx F(t>a/A(Xa,Q 2 )(t>b/B(xb,Q 2 ) X F cc^J /4,(q 2 ) , (1) 

a.b 

where ^ a , runs over all parton flavors, and Q 2 = q 2 + Am 2 . The K factor takes into account the higher order 
corrections. The incoming parton momentum fractions are fixed by kinematics and are x a — {y~Xp + 4Q 2 /s + xp)/2 
and Xb = {yj x'p + 4Q 2 /s — xp)/2. Quark annihilation and gluon fusion are the major sub processes for cc production. 
In the leading log, they are given by Jic| l, 

^ CE (Q 2 ) = §^(1 + 1)^/1^, (2) 



= gi(l + 2 + 1.) log(i±^) - ( 7 - + g 7)v ^ L (3) 

(4) 

where a s is the QCD running coupling constant and 7 = Am 2 /Q 2 . In Eq.[l] F ce ^j/^(q 2 ) is the transition probability 
that a cc pair with relative momentum square q 2 evolve into a physical J / tp meson. Qiu, Vary and Zhang considered 
three different parametric forms (representing different physical processes) for the transition probability. All the three 
forms could describe the experimental energy dependence of total J / tp cross section in hadronic collisions |0] . 

In a nucleon- nucleus /nucleus- nucleus collision, the produced cc pairs interact with nuclear medium before they exit. 
Observed anomalous nuclear enhancement of the momentum imbalance in dijet production led Qiu, Vary and Zhang 
H to argue that the interaction of a cc pair with nuclear environment, increases the square of the relative momentum 
between the cc pair. As a result some of the cc pairs might gain enough relative momentum squared q 2 to be pushed 
over the the threshold to become open charm mesons. Consequently, the cross sections for J/tp production are reduced 
in comparison with nucleon-nucleon collisions. If the J/tp meson travel a distance L, the transition probability F cS {q 2 ) 
in eq.|l| will be changed to 

Fcc^j/^q 2 ) -> F c ^j^(q 2 + e 2 L), (5) 

with e 2 being the square of relative momentum gained by the cc pair per unit length of nuclear medium. Of the three 
different parametric forms of the transition probability, all of which fitted the energy dependence of the J/tp cross 
section in hadron- nucleus collisions, only the following form, 

F^j/^q 2 ) = N JN ,e{q 2 )(l - 4m /_ 4m2 r*> (6) 

could describe the experimental J/ip data as a function of effective nuclear length |9|. For completeness purpose, we 
have redone the calculation of J/tp production as a function of effective nuclear length. We have used the CTEQ5 
parton distribution functions In Fig.l, NA50 data Q on J/tp cross section for proton- nucleon, proton- nucleus 
and nucleus-nucleus collisions as a function of the effective nuclear medium length L(A,B) is shown. The solid line 
is a fit obtained in the model. The parameter values, KNj/^=A58, e 2 =.225 GeV 2 / fm and of=1, are very close to 
the values obtained in Ref. |J. In the next section, we will use these parameters to analyze the NA50 data on the 
transverse energy dependence of J/tp to Drell-Yan ratio. 



III. E T DEPENDENCE OF J/tp IN PB+PB COLLISIONS 

NA50 collaboration presented transverse energy dependence of J/tp to Drell-Yan ratio in Pb+Pb collisions [Q. 
As mentioned in the beginning, the data shows anomalous suppression, which goes beyond the conventional nuclear 
suppression. In the present section, it will be shown that, the data are fully explained in the model of Qiu, Vary and 
Zhang, which treat J/tp suppression in nuclear environment in an unconventional manner. 

The Drell-Yan pairs do not suffer final state interactions and the cross section at an impact parameter b as a 
function of Et can be written as, 
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d 2 a DY /dE T d 2 b = o% J d 2 sT A (s)T B (s - b)P(b, E T ), (7) 

where a^Jj is the Drell-Yan cross section in NN collisions. All the nuclear information is contained in the nuclear 
thickness function, Ta, B (s)(= / (Izpa,b(s, z). Presently we have used the following parametric form for pa{t) ||, 

PA{r) = P ° f r-r , ( 8 ) 

with a = 0.53/m, tq = 1.1A 1 / 3 . The central density is obtained from J pA(r)d 3 r = A. In Eq.0, P(b,Ex) is the 
probability to obtain Ex at an impact parameter b. Geometric model has been quite successful in explaining the 
transverse energy as well as multiplicity distributions in AA collisions Jl^Jl^j . Transverse energy distribution in 
Pb+Pb collisions also could be described in this model B. In this model, Et distribution is written in terms of 
Et distribution in NN collisions. One also assume that the Gama distribution, with parameters a and j3 describe 
the Et distributions in NN collisions. Pb+Pb data on Et distribution could be fitted with a — 3.46 ± 0.19 and 
(3 = 0.379 ±0.021 @. 

While Drell-Yan pairs do not suffer interactions with nuclear matter, the J /ip mesons do. In the model of Qiu, Vary 
and Zhang || , suppression factor depend on the length traversed by the cc mesons in nuclear medium. Consequently, 
we write the J/ij) cross section at an impact parameter b as, 

d 2 a J/ ^/dE T d 2 b = a J J% J d 2 sT A (s)T B (s - b)S(L(b, s))P(6, E T ), (9) 

where er^y is the J/ip cross section in NN collisions and 5(L(b, s)) is the suppression factor due to passage through 
a length L in nuclear environment. At an impact parameter b and at point s, the transverse density can be calculated 

as, 

n(b,s) = T A (s)[l - e - CT ««T B (b- s )] + Ts ( b _ s )[! _ e -«7 JVW rAW] ) ( 10 ) 
and the length L(b, s) that the J/ip meson will traverse can be obtained as, 

L(b,s) = n(b ) s)/2p (H) 

Suppression factor S(L(b, s) can be calculated using Eq.^, noting that cc pairs gain e 2 momentum per unit length 
L. Parametric value of e 2 , as shown before was obtained by fitting nucleon- nucleus and nucleus- nucleus J/tp cross 
section data containing all Et- However, Eq.^| corresponds to a particular Et- Accordingly, momentum gain factor 
e 2 , needs to be modified, We modify the momentum gain factor e 2 to take into account the Et dependence as, 

where Eq is the momentum gain factor for all Et (which was obtained by fitting experimental data). L(Et) is the 
length through which a J/tp meson with transverse energy Et will travel. The length L(Et) can be calculated 

fd 2 bd 2 sTA(s)T B (b-s)[TA(s)+T B (b-s)]P(b,E T ) 
[ T) 2a NN p Jd 2 bd 2 S T A (s)T B (b-s)P(b,E T ) { 1 

Fluctuations of transverse energy at a fixed impact parameter plays an important role in the explanation of the 
NA50 data. Above 100 GeV, i.e., approximately at the position of the knee, the 2nd drop in the data is due to the 
fluctuations in Et ■ In order to account for the fluctuations, following Capella et al [Q, we calculate, 

F(E T ) = Et/E^ f {E t ), (14) 

where, 

f nf (f _ Jd 2 bEr(b)P(b,E T ) 

Et {Et) ~ Jd%P(b,E T ) (15) 

The function F(Et) is unity up to the knee of the distribution, and increases thereafter, precisely, where fluctuations 
dominates. The replacement, 
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L(b,s) ^ L(b,s)F(E T ), 



(16) 



then properly accounts for the fluctuations in the Et distributions. 

In Fig. 2, we have compared the Et distribution of J/tp to Drell-Yan ratio, obtained in the model with the exper- 
imental data obtained by NA50 collaboration. In the calculation, we have used <jnn =32 mb and /&nn —53.5 
We obtain excellent agreement with data. The second drop at £+=100 GeV is correctly reproduced. It may be 
emphasized that the present calculation is essentially a parameter free calculation. The few parameters of the model 
were obtained previously from the fitting the total J/tp cross section in pA and AA collisions. Excellent agreement 
with data indicate that the NA50 data is fully explained in terms of suppression in nuclear environment. 



IV. PREDICTION FOR RHIC ENERGY 



Present model can be used to predict Et dependence of J/tp to Drell-Yan ratio at RHIC energy. Recent PHOBOS 
experiment Jl5| showed that for central collisions, total multiplicity is larger by 70% at RHIC than at SPS. Et can 
be assumed to be increased by the same factor. Accordingly, scaled the Et distribution for Pb+Pb collisions at SPS 
energy can represent the experimental Et distribution at RHIC energy for Au+Au collisions (small mass difference 
between Au and Pb can be neglected). We have fitted the rescaled Et distribution in the geometric model to obtain 
the parameters, a=3.09 and /3=0.495 p6| . Nucleon-nucleon inelastic cross section ((Inn) was assumed to be 41 mb 
at RHIC, instead of 32 mb at SPS @. 

At RHIC energy the so-called hard component which is proportional to number of binary collisions appear. Model 
dependent calculations indicate that the hard component grows from 22% to 37% as the energy changes from y / s=56 
GeV to 130 GeV |J]. J/tp suppression will strongly depend on the hard component, as it effecti vely increases the 
density of the nuclear medium. For / fraction of hard scattering, transverse density n(b,s) in EqJlOj is modified to 

n mod {b, s) -» (1 - f)n(b, s) + fn hard (b, s), (17) 

with n hard (b, s) = <7nnTa(s)Tb (b — s). With hard component, transverse density is increased, as a result, suppression 
will be increased at RHIC energy. In Fig. 2, the thick solid line is the prediction for J/tp to Drell-Yan ratio at RHIC 
energy, for Au+Au collisions, obtained with 37% hard scattering component in the density. Very large suppression 
is obtained. Effect of Et fluctuations is not visible anymore (very large suppression washes out Et fluctuations). 
It is interesting to compare the present prediction with other model calculations. In fig. 2, the thin dotted line is 
the prediction obtained by Blaizot et al. [jl7| in model where all the J/tp mesons melts above a threshold density, 
essentially in a deconfined scenario. Very close agreement between the predictions obtained in a nuclear environment 
and in a deconfined scenario is interesting. It seems that it may not be possible to confirm the deconfincmcnt phase 
transition, which is expected to occur at RHIC energy, from the J/tp data. Recently several authors have proposed 
that at RHIC energy, in a deconfined scenario, recombination of cc pairs will lead to enhancement of J/tp's, rather 
than its suppression |l9f| . Inclusion of recombination effects may mask the large suppression obtained by Blaizot 
et al. fTi]] . However, nuclear suppression as calculated presently will remain unaltered. It may then be possible to 
distinguish the deconfinement phase transition from the J/tp data. 



V. SUMMARY 



To summarize, we have analyzed the NA50 data on transverse energy distribution of J/tp to Drell-Yan ratio in 
Pb+Pb collisions. The data were analyzed in a model, where suppression of J/tp is due to gain in relative square 
momentum of cc pairs as it travels through the nuclear environment. Some of the cc pairs can gain enough momentum 
to cross the threshold to become open charm mesons. The model, without any free parameters can well explain the 
NA50 data on J/tp suppression. Present analysis clearly shows that it is not essential to assume a deconfined scenario 
to explain the NA50 data. The model was used to predict Et distribution of J/tp to Drell-Yan ratio at RHIC energy. 
At RHIC hard component of scattering may be important. Very large suppression is obtained if the hard component 
is included. Interestingly, the prediction obtained in the model, with only nuclear suppression agrees closely with 
the prediction obtained in a deconfined scenario. However, as suggested by several authors, recombination of cc in a 
deconfinement scenario may lead to enhancement, rather than suppression of J/tp at RHIC energy. Recombination 
effect will not affect the nuclear suppression. Observation of enhanced production J/tp at RHIC energy will then 
confirm deconfinement phase transition. 
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FIG. 1. Total J/ip cross sections with the branching ratio to /i H 
collisions, as a function of the effective nuclear length L(A,B). 
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FIG. 2. Open and closed circles are the J/ip to Drell-Yan ratio in a Pb+Pb collision obtained by NA50 collaboration in 1996 
and 1998 respectively. The thin line is a fit to the data in the model described in the text. The thick solid line is the prediction 
obtained for Au+Au collisions at RHIC energy, with 37% hard scattering component (see text). The thin dotted line is the 
prediction obtained by Blaizot et al. |nj for Au+Au collisions at RHIC energy, in a model where all the J/ip mesons melts 
above a threshold density. 



